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Institute for Materials Science, 2-1, Sengen 1-chome, Tsukuba*shi, 
Ibaraki 305 0047 Japan, one of the inventors of USSN. 09/926,188 do 
hereby solemnly and sincerely declare : 

(1) THAT the attached is the REFERENCE FIGURES AND 
ITS EXPLANATION to prove the substrate surface can be obtained in 
which steps of each mono atomic layer are arranged in series by, in 
accordance with the present invention, polishing a diamond (100) 
substrate in the range of 1.5 to 6 degrees, and hydrogen plasma 
treating said substrate. And 



explanation that the diamond thin film growth takes place by 
step-flow growth mechanism by using this substrate, and therefore 
the S dope n type diamond of the present invention has the crystal 
perfect-ness which can not be realized by prior art. And 

(3) THAT the above-mentioned REFERENCE FIGURES are 
the graphics made from the data measured by AFM(inter-atomic force 
microscope) treating on Sep 10, 2003 by myself by the described 
conditions in order to prove that the substrate of the present 
invention possesses the crystal perfect-ness which can not be realized 
by prior art. 

And, I, Toshihiro Ando, being sworn state that the facts set 
forth above are true. 



(2) THAT this attached document is for the purpose of the 
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Molecular-orbital theory of monatoraic and diatomic substitutional defects 
as shallow n-type dopants in diamond 
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A survey of potential shallow n-typc donon in diamond was performed using dm atom superposition and 
electron ^localization molecalaiworbitel method. Electronic structure and defect stabilities were estimated for 
N, F, O, S, F, and Q in substitutional monovacancy sites; I, Xe, and the atom pairs BO. BS, NO. NS, and Bed 
were studied on substitutional divacancy sites. The calculations Indicate that the BS, NS, and Bed pairs aTid 
I in divacancy sites may provide shallow donor levels, and may be mora stable than P in a monovacancy site, 
[S0163-1 829(96)05 1 44*2] 



I* INTRODUCTION 

The experimental search for shallow «.rype diamond con- 
tinues. It is well understood that a main group atom with five 
valence electrons will, when substituting for a carbon atom 
in the diamond lattice, be a donor. 1 * 2 Substitutional nitrogen 
i$ known from measurement to be a deep optical donor (ab- 
sorption edge at 2 eV) and a deep thermal donor (activation 
energy 1.6-1.7 eV). 3 Theoretical calculations concur, 1,4 ' 3 and 
Show how one C-N bond increases in length from the dia- 
mond value of 1-54 A to about 1.99 A, while the other three 
bonds are compressed to about .1.48 A, in a trigonal distor- 
tion. Phosphorous, the neighbor Of nitrogen in the period 
below, is predicted to be a shallow donor, 4,5 This is due to 
the large size of the 3/> valence arbitals mat overlap the 2p 
orbitais of the neighboring carbon atoms strongly, even In 
the relaxed structure. The P-C antibonding &* orbital that i$ 
occupied with the extra electron can, at most, drop only a 
little beneath the bottom of the conduction band, and the 
trigonal distortion is small (in Ref. 5 three P-C bond dis- 
tances are calculated to be increased by 0.06 A, and the 
fourth by 0.09 A). This makes phosphorous potentially inter- 
esting as a substitutional n-type dopant in diamond, but it 
will be difficult to incorporate it because it id relatively un- 
stable. In Rcf, 5 ft was calculated to bind in a diamond 
monovacancy site about 10 eV more weakly than a carbon 
atom, whereas a nitrogen atom was calculated to bind only 4 
eV more weakly man the carbon atom. The weakness of the 
P-C bonds is caused by the promotion of the phosphorous 3p 
electron from its level, which is close to the bottom of me 
diamond band-gap region, to the top of the gap region. An . 
additional factor is the energy of compression of surrounding 
P-C and C-C bonds from their single bond equilibrium val- 
ues. Because of the expected relative instability of phosphor- 
ous in diamond, it is not surprising that attempts to incorpo- 
rate it by introducing compounds into the low-pressure 



growth medium have not yielded encouraging results."** 8 Re- 
cently, however, it has been reported that mild cold implan- 
tation of phosphorous ions followed' by rapid annealing 
. yielded n-type diamond with a low 0.2-eV activation energy 
for electron promotion to the conduction band. 9 This gives 
some hope for phosphorous. Also interesting -Is the recent 
report that phosphorous will codope with nitrogen during 
low-pressure diamond growth. 10 Such systems have, how- 
ever, not demonstrated n-type conduction or phosphorous- 
related luminescence In these systems the nitrogen concen- 
tration is about equal to or greater than that of phosphorous. 
If both dopants exist as substitutional atoms, perhaps all sub- 
stitutional P are oxidized to P + . An equal number of substi- 
tutional N would be reduced to N" with two band-gap elec- 
trons to give a longer C-N bond &* orbital with an energy 
even deeper in the band gap. Alternatively* substitutional PN 
dimers f if they form, might be deep donors. A theoretical 
study using the methods of this paper has, in fact, supported 
this interpretation of the lack of electrical activity and 
luminescence, 11 

Arsenic has also been implanted in diamond' 2,13 and Ru- 
therford backscattering spectroscopy indicated that —40- 
, 50 % of the As atoms were located in substitutional sites. 
Such structures should be highly stressed because of the 
large size of die As atom, and they will be donor defects. 
Measurements of the electrical behavior indicated the forma- 
tion of a semiconducting surface layer with a 0.41-eV acti- 
vation energy., 14 However, similar results were obtained for 
the implantation of Ne and Na, Indicating that lattice defects 
caused by the implantation may be responsible for the ob- 
served electrical properties. . 

Interstitial atoms, except for H; which has been predicted 
to act both as a deep donor and a deep acceptor, are ex- 
pected to be shallow donors in diamond The large overlap of ' 
the interstitial atom orbitais and the neighboring carbon atom 
orbitais result in substantial qestabilization of the occupied 
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TABLE I. Atomic parameters used Id the calculation* Principal 
quantum numbers /i, Slater orbital exponents f (a.u.), and diagonal 
Hamlhomsn matrix elements H (eV). 



Atom 


n 


( 




n 


f 


ff 


H 


1 


1.2 


-13.6 








Be 


2 


1.0050 


-10.822 


2 


1.0060 


-8.100 


B 


2 


1.2881 


-13.930 


2 


1.2107 


-9.298 


C 


2 


U174 


-16.590 


2 


L7717 


-11.260 


N 


2 


1,9237 


-18330 


2 


1.9170 


-11530 


0 


2 


12458 


-29.480 


2 


12266 


-14.620 


F 


2 


15635 


-39.850 


2 


2J500 


-19.420 


Al 


3 


1.3724 


-10,620 


3 


L3555 


-5.986 


P 


3 


1.8806 


-17.650 


3 


1.6288 


-11.990 


S 


3 


2.2223 


-22.200 


3 


1.8273 


-11360 


a 


3 


13562 


-26.540 


3 


243388 


-15.010 


i 


5 




-24,«!0 


5 


16790 


-14.450 • 


Xe 


5 


3.1433 


-23.700 


5 


18439 


-13.440 




cr* orbitals, which will lie near or above the bottom edge of 
the conduction band, despite the repulsive forces chat serve 
to push the neighboring carbon atoms off their lattice sites. 4 * 3 
Thus quantum calculations predict that interstitial Li (Reft. 4 
and 5) and Na (Ref. 4) are shallow donors. By contrast, 
. substitutional Lf is predicted to be a shallow acceptor * Pho- 
toemissioa measurements with indiffused Li show the pres- 
ence of acceptor levels 1.0-1-5 ©V above the valence-band 
edge, 6 These levels might be due to Interstitial H/ 3 substitu- 
tional Ll 3ao or something else, but currently there is no 
proven answer, and experimental investigations continue, 2 

IL RATIONALE FOR SYSTEMS STUDIED HERE 

We examine, the properties of selected atoms in monova- 
cancy and divacancy substitutional sites and selected pairs of 
atoms in divacancy sites, with a view toward identifying 
shallow donors for diamond Defect .stabilities, electronic 
structures, and the effects of structural relaxation on optical 
and thermal excitation energies are calculated. For the mono- 
vacancy substitutional sites we consider N, P* O, S, F, and 
G, Oxygen and sulfur, from group VI B, have two electrons 
that will be promoted into the band gap, and F and Ci f from 
group YIIB t have three. Whether these are deep or shallow 
donors should depend on the electron count and the orbital 
sizes, 1 * 9 and this study addresses these issues. Although they 
were studied previously, N and P, from group V B, are in* 
eluded for comparison. 

We also studied the incorporation of large single atoms 
and pairs of atoms into divacancy sites. A divacancy site has- 
six dangling orbitals extending into the cavity, compared to 
four for a single-vacancy site, and there is one electron for 
each orbital. Such vacancies could act as traps for interstitial 
H» as we discussed in a previous paper. 15 In the present con- 
text we are interested in characterizing the binding of atoms 
within these divacancy sites. We consider the resulting de- 
fect electronic structures In terms of orbital correlations of 
the dangling vacancy orbitals with the orbitals of atoms in- 
troduced into the vacancies, as In our earlier work involving 




FIG. 1, Bulk superimposable C71 and cluster models used 
for defect studU*. The cross hatched central atoms are removed to 
create vacancy sites or replaced by foreign atoms to create substi- 
tutional defects. Neares>neighbor (shaded) atoms around the de- 
fects are allowed to relax in die (111) directions. The dangling 
lurfaco^wbitaJ-termlnadng H atoms aw not shown for the C« clus- 
ter. 



substitutional B and N. 1 Just as an atom with five valence 
electrons, such as N and P wih\ when in a monovacancy site* 
promote one electron into the band gap. an atom with seven 
valence electrons, such as I, will promote an electron into the 
band gap when in a divacancy site. Whether such an I atom 
is a shallow donor depends on the strengths of the orbital 
overlaps. We have explored thia issue for I and also Xe. 
which will promote two electrons into the gap when in a 
divacancy site, Finally, for a diatom pair placed in a diva- 
cancy site, the electronic structure can be thought about in 
two ways. For example. BO is isoeleetronic to CN, and so 
adding it to a divacancy site is similar to adding N to a 
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FIG. 2, Structure variable* for the monovacancy and rnonoxnb- 
. arftutlonal mom site. 

monovacancy site. Alternatively, if we recognize that, when 
present in the diamond lattice, there will be a single c bond 
between B and Q t the remaining seven valence electrons will 
be in orbitals that will form bonding and aotibonding coun- 
terparts with the six diYocancy dangling orbitals containing 
six electrons, thus promoting an electron into the gap. In 
either viewpoint, one electron will be promoted into the band 
gap or above, and the specific properties of the atom will 
determine whether the defect is a shallow donor, la this re^ 
gard we examine BO, B$ ( and BeCl as potential donors with 
one electron in a single orbital promoted into the band gap, 
and NO and NS with three electrons in two orbitals pro- 
moted into the gap. 



Ttt THEORETICAL APPROACH 

The atom superposition and electron deJocalizadon mo- 
lecular orbital (ASED-MO) theory 17 was employed for this 
study. This is a semlemplrical theory for calculating struc- 
tures and electronic properties of molecules and solids from 
atomic data, namely, valence atomic orbitals and Ionization 
potentials. The theory is based on partitioning the electronic 
charge density distribution function into atomic and deiocal- 
i*ation (doe to forming bonds) components p fl and p rfl re- 
spectively, 



(1) 




- stretched 



PIG. 3. Shapes of occupied or* orbitals for trigonalJy distorted 
N, 0, and S substitutional defect centers (o, orbital on left), and for 
tctragonally distorted F and CI substitutional defect center? (b 7 and 
a, orbitals Ott right). 

Pairwise atom superposition energies E* b are calculated by 
integrating the electrostatic force on the nucleus of the less 
electronegative atom of each pair of atoms a and b t as they 
are brought into the molecular configuration. The sum E R is 
called the atom superposition energy, 



a<b 



(2) 



The electron defooalittition energy E D is the sum of integrals 
of the electrostatic force on the same nuclei as one atom after 
another is added to the cluster. This cannot be evaluated 
because the functional forms of p d are generally not avail- 
able, but E D is sufficiently well approximated by a 
molecular-orbital electron derealization energy Ae mo that' 
the total energy 



can be replaced by 



E~B R +AE MO 



(3) 



(4) 



for many applications. The AE MQ that works in Eq. (4) is 
calculated using a modified extended HucJkel HaoiUtonian 
where diagonal matrix elements are set equal to the negative 
of measured atomic valence state ionization potentials, 
sometimes shifted for highly ionic systems. Off-diagonal el- 
ements for orbitals on the Same atom are set to zero, and for 
orbital i on atom a and orbital j on atom b they are given as 

^lA^{^Hp^~QA2R ab )S^ (5) 



TABLE IL Calculated Wndtog energies, Bfi (eV), for atoms to vacancy sites and structure parameters as 
defined In Fig. 2 for the optimized trigonal and tetragonal defect strnctures; All displacements (z and d) and 
interauclear distance* (R) arc in A units. Atoms that are next-nearest neighbor to substitutional atoms A are 
constrained to the bulk diamond lattice sices. 



Trigonal 



Tetragonal 





Atom 


BE 












. BE 




dtx 


du 


&A-iJl 


Ra-x* 




N 


10.7 


0O33 


0.224 


-0.007 


\S99 


I.47I 


9.60 


0.189 


-0.006 


0.135 


1.791 


1.434- 




0 


7.05 


0.347 


0.216 


-aoi i 


2.113 


1.450 


6,09 


0355 


0.234 


0.007 


2.000 


1373 


1 


F - 


2.67 


0.2*57 


0.236 


0326 


2044 


1.795 


3.66 


0.507 


0.214 


0.075 


Z088 


1386 




P 


4.78 


O004 


aa&o 


0.063 


1.624 


1.602 


4J8 


0.003 


0.074 


0.061 


1.600 


1.615 




S 


-5,67 


ai§4 


0.475 


0,186 


1199 


t.674 


-734 


0.020 


ano 


0.141 


1.722 


1.670 




a 


-10,03 


0.420 


0.223 


0.144 


2.183 


1.623 


-9.44 


0.255 


0.425 


0.189 


zm 


1.595 
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FIG, 4. Bond-gap electronic structures calculated for the va- 
cancy and for monosubttitmkmal atoms in the clusters, The 
filled valence band la shaded and the empty conduction band Is 
unshaded. Occupied band-gap orbital energy levels are shown for 
the vacancy V. which has four dangling orbltals with four electrons, 
and for the dopant atoms N through CL For N, O, F, S, and CI 
orbital* that are bonding between valence s and the cluster orbitals 
arc shown beneath the valence band in energy. 

where R wh is the mtetrtucLear distance, 5?/ the overlap inte- 
gral, and the /Ts on the right-hand side are the diagonal 
matrix dements defined above. 

The single-f Slater orbital exponent and diagonal matrix 
elements used in this work arc(bTab^DThe carbon atom 
orbital exponents arc those cboseTrtcTpfScluce the C-C bond 
distance in graphite {1.42 A) using the ASED-MO band 
program. 1 * The carbon atom diagonal If matrix elements are 
the negative of the measured 2s and 2p atom ionization 
potentials. With these parameters the calculated C-C bond 
distance from ASED-band calculations for diamond was 
1.53 A compared to the measured 1.54 A. In the cluster 
models of this paper 1.54 A is used for C-C distances except 
when there is relaxation included around the defects being 
studied. These carbon atom parameters yielded 8.05 eV per 
atom for the atamization energy of diamond in the ASED* 
band calculations* compared to an experimental determina- 
tion of 735-eV per atom. The bulk modulus was calculated 
to be 539 GPa, which was within the range of experimental 
estimates, 442, 556, and 575 GPa, These results were similar 

i rejaxation £rvbx> and 
(thermal E Ul _jJtetivatioti energies pSyfiar prornodnglar^ ectron 
to the eonduSBon Bnct**'" "" rr ~ " ' ■ ■ « »■ - 



Defect 



N 


3.57 


2.63 


0.94 


0 


4,74 


1.22 


3.51 


F 


4.90 


1.19 


17! 


P 


-<UI 


1.14 




S 




2.07 




Cl 


2.97 


3.71 


-0.79 ~ 


ao 


4.11 


1.18. 


2.92 


BS 


-0.12 


0.81 


-C.93 


NO 


4,04 


1.64 


238 


N5 




3.12 


-0.03 


Bed 


X41 


2.31 


0.1 i 


I 


0.95 


1.09 


-0.14 


Xe 


452 


ass 


3.67 



TABLE rv. Calculated stiucuirc parameters rot A %). See Fig, 
2 and caption to Table IL 



Atom 




d, 








N 


0.000 


-0.150 


-0.015 


1.525 


1J25 


O 


0.277 


0.242 


-0.016 


2.0SS 


1.456 


F 


0314 


0.216 


0.265 


2.070 


1.726 


P 


aooo 


0.070 


0.070 


1.610 


1.610 


S 


-aoas 


0.147 


0.1*5 


1.683 ' 


1.697 


a 


0229 


0419' 


0.141 


2.188 


J.619 



to those calculated by density-functional band theory, IM 
A, 7.94 eV. and 43S GPa for the bond length, atoraization 
energy, and bulk modulus, respectively. 19 Generally, the 
single-f orbital exponents are chosen to flt molecular di- 
atomic bond lengths. Thus the theory is parametrized for 
accuracy for bonded structures. The parameters for Q N, and 
P were used In Ref. 5, which can be consulted for further 
details. 



rV. DIAMOND MODELS 

We have used a V*CtoH m cluster to model the single va- 
cancy and an AC^H^j cluster to model a substitutional atom 
defect For the di vacancy and disubstiturJonal defect the 
V^CpeH?* and ASC^tl n cluster? (Fig. ]) were used. Here v 
represents a single vacancy and A and B represent substitu- '- 
tlonal atoms. All defect structures were optimized by allow- 
ing neighboring C atoms to relax toward or away from the 
center of the vacancy site along (111) directions. For the 
clusters used, no surface atoms are relaxed for they we next- 
nearest neighbors to the defects. Comparison with previous 
calculations, wherein next-neighbor relaxations along (ill) 
directions were allowed, showed 0.1-0.2 eV reductions in 
calculated binding energies of atoms to monovacancy sites, 
and in such cases cluster surface atoms were allowed to 
move. 3 We feel the present model, which is computationally 
more efficient, is comparable to or better than the previous 
one; A few calculations of relaxations allowing deviations 
from the (J 11) directions of motions of relaxing atoms re- 
sulted in small — 0.1-eV stabilizations that can be regarded 
as negligible at this level of modeling. Relaxation of atoms 
surrounding the monovacancy gives rise to a tetragonal dis- 




FIG. 5. Structure variables for divacancy and dUubfidtudonai 
sites. 
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FIG. 6. Band-gap electronic structure* for fiMeronuclear occu- 
pations of the adjacent divacancy site in the cluster and also for 
1 and Xe in this aita, See caption to fig, 4, the divacancy ( V 2 ) has 
six dangling orbitals and six electron!. 

tortion of C^ symmetry with two C atoms relaxed away, by 
0.187 A, the other two by 0.096 A. and a stability gain of 
0.36 eV. For the divacancy a stability of 0.89 eV comes with 
relaxation of the six surrounding C atoms to a structure with 
Oi h symmetry. This structure has an inversion center with* at 
each end, two C atoms relaxed away 0,094 A, and the third 
by 0.17 1 A. When the defect atoms are placed in the diva- 
cancy, the surrounding C atoms can experience large shifts 
in position: these shifts are related to the atom size and num- 
ber of valence electrons as discussed below. 

Energies for promoting band-gap electrons into the con- 
duction band were determined in two ways. The optical en- 
ergy was taken as the difference between the highest occu- 
pied band-gap orbital energy and the lowest empty 
conduction-band-orbital energy. The thermal energy was de- 
termined as the difference between the toial energy of the 
defect cluster, with the structures relaxed aa described above, 
and the total energy of the relaxed positively charged cluster 
plus the binding energy of an electron in the C-C bond 
stretched state as discussed in Ref. 5. There it is shown that 
the calculations favor localization of the promoted electron 
in this cluster and delocalization into the conduction band in 
the bulk. For cases of substitutional dinners the cluster is 
used for this, and for monomers the C? 1 cluster is used with 
a stretched bond involving displacement of the central atom 
and a neighbor so that the defect has the same local struc- 
tures as determined in Ref. 5 for the cluster. To three 
figures after tfte decimal point the calculated electron affini- 
ties of these two clusters are the same 

We also determined the bonding energy of A and B in the 
di vacancy site and isolated A and B in roonovacancy sites 



AE 

AB(s)-+AU)+8(s) ( 6 ) 

where s stands for substitutional. The energy change for Eq. 
(6) was calculated using the total energies E of C^H™, 
AC,^ SCtoH^ C»sH 7S , and ABC^i^. 

A^ = «Cpt) + fiCAW)+£0C*))-2£(C,,)--^aW>. 

(7) 

V, RESULTS AND DISCUSSION 
A. Substlra«Ottal atoms? N, P, 0, S> F, CI 

Though N and P were analyzed In the earlier paper/ the 
calculations were repeated using the slightly different con- 
straints of the present work for comparison and to establish 
trends. Structure parameters defining the substitutional atom 
defect structures are shown in Fig. 2. Tho central (substitu- 
tional) atom/i has been allowed to move up and down 
wi«3 respect to the lattice point at z=0, giving rise to 
trigonal distortions and along the lino in the xz plane (rf^) to 
give tetragonal distortion*. The four shaded surrounding at- 
oms were allowed to relax inward (-d) and outward (+d) 
in the <lil> directions. 

From the calculated structure and energy parameters in 
Table U, ft Is seen that only the atoms with three extra elec- 
trons, F and CI, favor the tetragonal structure, over trigonal. 
Atoms with two extra electrons, 0 and S, favor the trigonal 
structure. Of the atoms with one extra electron, N favors the 
trigonal structure and P shows no preference. . 

Focusing on the trigonal structure of N, O, and F, we sec 
for the first two a large elongation of one bond to 2,0 and 2. 1 
A and large compressions of the other three bonds to L41 
and L48 A (the diamond value is 1.54 A). The elongations 
are the consequence of the nearest-neighbor antibonding na- 
ture of orbitals in the conduction band. A bond stretches to. 
stabilize one of the N-C or O-C a* orbitals, but this occurs at 
the expense of Coinpressing three N-C or O-C bonds. The 
occupied 0* orbital is shown in Fig. 3 based on the 0 caU 
CuJaa'on. When a third extra electron is brought in with F, the 
trigonal structure is mfforaln ed. and the other three bonds 
stretch as shown CfaT Table jb However, more stability is 
achieved in the tetragonal structure with two Yery long bonds 
of 2.1 A and two short ones at 1.4 A. This allows two a* 
orbitals to drop deep into the band gap, one for a doubly 



TABLE V. Calculated bonding energies for atoms A and B in the relaxed divacancy site, BE(V 2 ) (eV), and In monovacancy sites 
BE(2 V) (eV), ABU) dhwociarJon energies A£ (eV) as defined In Eqj. (6) and (7), and structure parameters as defined in Fig. 4 for the 
optimized defect structures. All displacements from lattice points d, and internuclear and otber distances R are in A and angles in degrees. 
Carbon atoms that are nexi-ncarest neighbors to substitutional atoms A and £ are constrained to the bulk diamond lattice sites. 



Dlmer 


Bfi(V 3 ) 








a 




0 


*a 






<*4 










*a-3,« 


BO 


17.92 


18.28 


3.09 


1482 


9.8S 
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TABLE VL Calculated structure parameters for A See Fig. 4 and caption to Table V. 



Dimer Ram a R* 0 .R$ d x d-# rf 4 *w A/h *j-a #fl-jjtj 

BO 4 1.662 0.00 1340 0,00 1.602 -0.001 -0.001 0.000 0.000 1,539 1.539 1.521 1.521 

BS + 1.60S 050 1.678 0.50 1.526 0.040 0.044 0,156 0.165 1.554 1536 1.713 1.704 

NO* 2.114 0.00 U62 0.00 1.791 -0.004 -0.004 -0.015 -0.016 1.472 1.472 1,450 1.450 

NS + 1.507 0.03 U*4 12.30 1.451 0,275 O.W 0.156 0.136 1134 1.454 1.729 1.708 

BcCi* 1.608 0.00 1525 0,00 1.623 0,045 0.045 0.113 0.U3 1.559 1.559 1.658 1.65$ 



occupied and the other for a singly occupied orbital. These 
orbitals are shown In Fig. 3 based on the CI calculation; For 
Ci the C\ orbital is singly occupied, and the b 2 orbital is 
doubly occupied; for P the order Is reversed. 

Moving now to P, S, and CU ihe near equality of stabili- 
ties tor the trigonal and tetragonal structures of substitutional 
P reflects the large size of the P atom. The lengthening of the 
P-C bond for the trigonal structure is small compared to our 
finding for substitutional N, and the other three bonds are 
shortened only slightly. Being smaller, S does not have such 
a strong isotropic component to its force against the lattice 
(prior to relaxation) and for the trigonal Jtructure one bond is 
able to stretch considerably to 2.2 A* accommodating the two 
band-gap electrons in a S-C cr* orbital, while die other three 
S-C bonds compress to 1.45 A, Substitutional CI favors the 
tetragonal structures with two bonds quite long at 2, 1 A, and 
the other two at 1.6 A As Rg, 4 shows, the sequence of 
donor levels lies higher for P, S, and CI than for N, O, and F, 
and this is dominated by the larger orbital sizes of the 
former, though reduced atomic ionization potentials, which 
affect the ASED-MO parameterization, also contribute. 

The optical excitation enex^esrbase^ on the energy levels 
shown in Fig. 4, are listexf^oj^e^wilong with the relax* 
ation energies based on opdntm^ffie AC 7 o + clusters and 
the C9»~ cluster. and the thermal excitation energies, which 
are calculated as the sum of the optical and relaxation ener- 
gies. Of the set, only P is predicted to be a shallow optical 
donor, and both P and CI are shallow thermal donors. 
The calculated relaxation energies are relatively large for 
S. and Chan d in fact^forQU^ r elaxation eT tergyj3U*se& 

>>w a2 £nangrfror a^leeao^^ *° gES jow thermal dono rJhe 
reJaxaQblTel^^ orbital 

^-TBTSQ^^ away 
^Tj ^e^a^ond lattice jg &cture. Relaxed structures of the 

?^C?o ^usiet^ w ^ > . f he kjtial sta te 

SEru^urSTn ^leltS ft qw how large stability flauwLfl flBUiCr 

_ for mos eca sea wrteTeTnere are Faxge ^^ctjrr e cranges. ^or 
N T removaIjf^ 
"bon JponH e 

r5TATv^ffN'*^]nuig m "the tetrahedral site. This relaxation 
gives 2.63 eV of greater stability. Because P is larger than N, 
relaxations with substitutional P are smaller. For O, relax- 
ations of structure and energy are small because after ioniza- 
tion one electron still remaina in a O-C o* oibital, so the 
stretched O-C bond remains long. For S the two band-gap 
..electrons settle in a very stretched (2.20 A) S-Co* bonding 
orbital, and upon ionization the structure relaxes to one simi- 
lar to P with relatively small S-C bondlengm variations. 
When an electron is removed from F, which has two 2.09-A 
bonds in the tetragonal structure, a trigonal structure with a 
. Jong F-C bond forms. Substitutional CI experiences similar 



structure change when ionised, but with a large relaxation 
stabilization because of its large compressive energy against 
the surrounding lattice when two Cl-C bonds are long. 

B. Substitutional heteroatom pairs m di vacancy sites: 
BO, B$, Bed, NO.andNS 

Structure parameters for these systems, as defined in Fig. 
5, are given in Table V along with stabilities relative to iso- 
lated substitutional atom states. Calculated Optical excitation 
energies, relaxation energies, and thermal excitation energies 
are given in Table HL All five substitutional dimers are cal- 
culated to be stable with respect to dissociation [Eq. (6)]. 
That is, the energy changes, AE in Table V, are all positive. 
Binding energies of A+B in divacancy sites compared to 
isolated monoyacancy sites are higher for BS, NS, and BeCl 
and lower for BC and NO. The difference between these 
binding energies is another measure of the relative stabilities 
of the associated and dissociated substitutional defects. The 
dissociation process of Eq. (6) favors the associated pairs in 
part because V 2 has one more C-C bond man 2 V and in oar 
calculations the strength of this bond is 3.45 eV. The impli- 
cation of these results is mat when the dimers form during 
growth may are Hxely to be robust. 

The above results can be analysed using what we learned 
about substitutional atom defects in Sec. V A. The substitu- 
tional BO defect is isoelectronic with the substitutional CN 
defect, and since B and C are adjacent In the Periodic Table 
and relatively close in size, as a« O and N, deep donor 
properties are predicted just as for N, and the O-C bond is 
elongated just as the N-C bond is. The ionized defect BO + 
takes trigonal symmetry (Table VI) about B and O. The BS 
defect behaves like substitutional CP, with which it is iso- 
electronic and shares a similar size. Because of the large 
sizes of sulfur, a band-gap orbital, involving a stretched B-C 
or S-C bond does not form, and when ionized the structure is 
not quite trigonal. The NO defect is similar to CF in that it is 
calculated to be a deep optical and thermal donor. Substitu- 
tional F takes a tetragonal structure with two very long (2.09 
A) and two very short (1.39 A) bonds to C N is in a roughly 
trigonal orientation, with three N-C bonds of the same (1.47 
A) length and the NO bond long at 2.10 A; O has one long 
(2.02 A) O-C bond and two short (U6 A) bonds. With this 
pseudotetragonal distortion about O, the doubly occupied or- 
bital has become about t eV more stable (Fig. 6) than the 
trigonal O level In Fig. 4. When ionized, the structures about 
N and O are trigonal. Substitutional NS, being similar to 
CCU is calculated CO be a shallow thermal donor, like CC1> 
and a deep optical donor. It has a pseudotetragonal structure 
with two long 2D2-A N-C bonds, When NS is ionized, the 
structure becomes pscudomgonal with one long N-C bond 
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due to the large size of S. The ionized NO defect, on the 
other hand, has a readily recognizable trigonal structure. Fi- 
nally, BeCL which is similar to CP, 1* unlike it, a deep 
optical donor but, since the long (2,05 A) CCi bond can 
relax 0.50 A on ionization, there is a large (Z30 eV) relax- 
ation energy that renders Bed a potential shallow thermal 
donor. The ionized defect takes on trigonal symmetry. 

C Large atoms la dfvacancy sites I and X* 
As explained in Sec. J, I should promote one, and Xe two 
electrons into band-gap orbitals when in divacancy sites. 
This is seen in Fig. 6, and from Table m we see that both are 
expected to be deep optical donors, and I. but not Xe, is 
predicted to be a shallow thermal donor. With respect to the 
free atom states, I binds with a stability 0 f 1.42 eV, and Xe is 
unstable by 4.49 eV. I fits well in this site, with the inversion 
center of the V 2 site retained. The relaxation of the two C 
atoms at each end is increased by 0.06 A to 0.154 A and for 
the third C atom the relaxation is decreased by 0.103 A to 
0.068 A. Nearest neighbor I-C bond lengths are four at 2,085 
A and two at 2,001 A. These are close to the methyl iodide 
value Of 2.13 A. Around 1+ all C-I bonds are 2.013 A long 
and all C relaxations away are 0,068 A. Xe has moved weU 
off center, so that it is L774 A from three neighboring C and 
2.515 A from the other three. Being close to three carbon 
atoms satisfies its smaller atomic radius. It forms bonds with 
these C atoms and three electrons are sent to the three dan- 
gling half-filled orbitals centered 2.515 A away, so that they 
become fully occupied These orbitals are shown in Fig. 7 
along with iodine's donor orbital for comparison. 

VI. CONCLUDING COMMENTS 

Our study provides a motivation for attempting to create 
shallow donors in diamond by incorporation of I and CI, and 
the pairs B with S, N with S, and Be with a It may be 
possible to form them by kinetic trapping during low- 
pressure growth or by ion implantation. 

Generally, we find that large substitutional atom size is a 
prerequisite tor shallow thermal doping. Half-filled donor 
levels are more likely to be shallow donors because they are 
destabilized more, being antlbondiog counterparts in situa- 
tions where there is a bond order of J. Double occupation 
means zero bond order, and the atoms can move farther 
apart, stabilizing the donor orbital until tbe restoring force of 
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FIO. 7. Donor orbital for I in the center of Iho divacancy site 
and a* and <f-»ymmetry band-gap orbitals for Xe in the off-center 
position within the dvaeaacy site. 

the lattice prevents further separation. Thus P, with one 
band-gap electron, Js predicted to be a shallow donor, and S, 
with two, is a deep donor, whereas the small atoms N and O 
are both deep donors. The mree»band.gap-electron atom CI 
is calculated to be a shallow donor because of a large relax- 
ation from tetragonal to trigonal structure upon ionization, F, 
on the other hand; is a deep donor because of its small size 
and relatively stable occupied band gap orbitals. 

For the substitutional dimers, results are similar to chose 
obtained for tho monomers, Each of these dimers has at least 
one first row atom. Thus BO is like N, has one band-gap 
electron and is a deep donor, while iscelectconic BS is iifce'P 
and is a shallow donor, as is also the case for BeCL The- 
three-band-gap-electroa dimer NO is like F and is a deep 
donor, whereas NS is, like CI, a shallow donor. 

. Finally, I in the dlvacancy site has a single electron high 
in the band gap and is a shallow donor. Its large size is 
responsible. The donor orbital containing two electrons for 
Xe is low in the gap because of this atom's small site and it 
is a deep donor* 



ACKNOWLEDGMENTS 

This work was supported by the National Science Foun- 
dation, Oram No. DMRr03527, and by the Chemistry De- 
partment of Case Western Reserve University. Discussions 
with Lubomir KostadinoY are acknowledged. 



*On leave from Institute of Physical Oiemistry, Bulgarian Acad- 

amy of Sciences, $ona, Bulgaria 1113. 
'S. P. Mehandru and A. B, Anderson, J. Mater. Res. 9, 3 83 (1994). 

This paper illustrates (he bonding principle and reviews the lit- 

eranue concerning substitutional N atoms. 
3 G. Popovid and M. A. Prclas. Diamond Kafcc Mater. 4. 1305 

(1995). ThU is a brief review of research toward synthesizing 

/i-type diamond from a pnenoznanological point of view. 
3 A. T. Collins and E C Ligmowlers, in Th* Pwpertiai qf Dia- 

mond, edited by J. & Field (Academic, London, 1979). 
4 S. A. Kajihara, A. Aatonelii, J. Buuholc, and R. Car, Fhyi. Rev. 

Lflit ^ 2010 (199i>. 
3 A. B. Anderson and S. p. Mehandru, Phys. Rev. B 48, 4423 

(1993). 



6 K. Okumara, J. Mart, and M. Machonkia. Appl. Phys. Latt 57, 
1907 (1990). 

7 1C Okano, JL Klyota, T. Iwas&ki Y. Nakamunu Y. Aldba, T. 
Karosa, M. gda, and T, Nakamura, AppL Phys. A 51, 344 
(1990). 

fl $. K Schauer. J. R. Flemish, R. WJttstruck, M. L Lands trass, and 

M, A, Piano, Appl. Phys, Lett 64. 1094 (1994). 
V. P, Prins, Diamond Reiat Mater. 4, 580 (1*995). 
l0 C. Z. Cao, F. A. J, M» Drtessea, Q. J. Banhufe, L. J, Oiling, and 

P. F. A. A&emada, J, Appl. Phys. 78, 3123 (19W), 
n A. B. Anderson and L. Kostadlaov, I Appl Phys, (to be pub- 
lished). 

"S. M Oorbatkia, R A. Zuhr, J. Roth, and H. Nararaoto, J, AppL 
Phya. 70, 2986 (1991). 



Received at: 4:54AM, 4/19/2004 

FROM ¥dj#&raJ?r/H I RAY AMA & CO- 0333522150 



2QQ4S 4fl1 9B <Ji) 1 8 : 1 9/TOl 8 : 0 4 /ggff-tS 3 0 0 9 0 3 8 4 7 P 22 



14 34* 



ANDERSON, GRANTSCHAROYA. AND ANGUS 



* Bhflmth.Ram, H, Quintel, M. Resile, C Ronning, H. Hofciisi, 
and 5. G. Jahn, J, AppL Phys. 78. 5180 (1995). 
J* D. Hunn, N. R. Parikh, M* L. Swansori, and R_ A. Zur, Dia- 
mond Relat. Mater. 2, 847 (1^3), 
15 S, P. Mehandni, A. B. Anderson, and h C Angus. J. Mater. Res. 
7. 689 (1992). 

■*tvf_ Rflstla, K. Bhftrum-Ram, H. Quintd. C. Ronning S. G. Jahn. 
and U. Wahl, AppJ. Phy*. Lait. 66, 2733 (1995). 



,7 A. B. Anderson, J. Che*,. Phy*. 62, J 1S7 (1975); A. B. Aadcrsoru 
R. W, Grimes, and S. Y. Hong, J. Phys. Chem. 91. 4245 (Z987); 
A* B. Anderson, Int J. Quantum Chein. 49, 581 (1994). 
K. Nath and A. B. Anderson, Solid State Commun. 66 277 
(1988). 

l9 K L Chan* and M. L. Cohen, Phys. Rev. B 35, 8196 (1987); M. 
T. Yin and M. U Cohen, ibid, 26, 5668 (1982). 



